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We examine the role of instantons in the zero-temperature chiral phase transi- 
tion in an SU(N) gauge theory. For a range of Nf (the number of fermion flavors) 
depending on N, the theory exhibits an infrared fixed point at coupling a*. As 
Nf decreases, a* increases, and it eventually exceeds a critical value sufficient to 
trigger chiral symmetry breaking. For the case N = 2, we estimate the critical 
values of Nf and a* due to instantons by numerically solving a gap equation with 
an instanton-generated kernel. We find instanton effects of strength comparable to 
that of gluon exchange. 
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Instanton configurations of the Yang-Mills potentials A a {x) have been studied 
extensively for over two decades. They play a central role in the solution of the QCD 
,iZ! U{1) problem [f2|, and a host of other physical consequences have been examined ||. 

In particular, many authors have studied their possible role in the dynamical breaking 
of chiral symmetry in QCD [|], ^, |5|, [7], Q. All these studies face a difficulty: their 
effects are dominated by large instantons, on the order of the inverse confinement scale 
of the theory, where the interactions become strong and instantons overlap. Reliable 
quantitative estimates are therefore difficult. 

A recent paper |J suggested that the chiral phase transition in an SU(N) theory at 
zero temperature, as a function of the number of fermions Nf, could be analyzed without 
the complications of confinement. For a certain range of Nf, the two- loop j3 function has 
an infrared stable fixed point, with the fixed point coupling a* increasing as Nf decreases. 
The transition was argued to set in when a* exceeded a certain critical value. That work 
considered forces arising solely from gluon exchange. 

Here, we examine the role of instantons in the same theory. The fixed point allows 
a more reliable study of instanton effects because it limits the growth of the effective 
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coupling at the large length scales which dominate the dynamics, better controlling the 
integration over the size of single instantons. We will begin our presentation with a brief 
review of the model; then display a gap equation with a kernel appropriate to an instanton 
background; next qualitatively consider the nature of possible solutions before displaying 
numerical results; and finally discuss reliability and draw some physical conclusions. 

We write the Lagrangian for SU(N) gauge theory as 

C = $(i$-g fl a T a )ip - jF^F atMJ + (gauge fixing terms 

where g is the gauge coupling, ip describes Nf flavors of Dirac fermions in the SU(N) 
fundamental representation, and the adjoint index a ranges over (l,---,iV 2 — 1). We 
assume a vanishing 9 parameter multiplying the anomaly term FF. 

After renormalization, the coupling a(fi) = g 2 (yu)/47r runs with energy scale, obeying 
a renormalization group equation 

d 
fjL—a(/j) = /3(a) = —ba 2 (n) — ca 3 (n) — . . . . (2) 

The first two coefficients are renormalization scheme independent: 

b = -L(11JV-2JV / ) (3) 
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UN Z - WNNf - 3 — T^N f . (4) 



Asymptotic freedom requires b > (or Nf < UN/ 2). If in addition c < 0, there is an 
infrared stable, non-trivial fixed point [|10|, located (to two-loop accuracy) at coupling 



b 

«* = --. (5) 

c 

Assuming that at some reference scale A, the running coupling has a value a (A) between 
and a*, we then have a(fi ^C A) — > a* and a(fi ^> A) — > 0. We can write the two- loop 
solution to the renormalization group Eq. (|2]) in the form 

-l/(6a») r -I 

A/yU = a*/a(/i) — 1 exp — l/6a(/x) , (6) 

giving a(A) = 0.7822 a*. Note that a(A) will vary when we change Nf/N. 

The reliability of this result depends on whether higher order terms in the (3 func- 
tion can be ignored. That is guaranteed for sufficiently small a*; more generally, the 
higher order terms are renormalization scheme dependent and so can be made arbitrarily 
small. For such special renormalization schemes, or for larger a*, there could of course 



be important higher order corrections to other quantities of physical interest Gj . Never- 



theless, the higher- loop j3 function JTT[] continues to display an infrared fixed point, with 
scheme-dependent accuracy discussed below. 

In Ref. H, it was argued that gluon exchange triggers dynamical chiral symmetry 
breaking when the fixed-point coupling a* exceeds a critical value 

ac 3C 2 (R) 3(iV 2 -l) ' l } 

This happens when Nf drops below the critical value for gluon exchange, 

The leading approximation for large N is quite good even for N = 2. 

An estimate of the higher order corrections to the gap equation driving the breaking 
provides some evidence |l^] for the reliability of these results. For a < a c , the next order 
corrections were found to be relatively small (less than 20%) compared to the ladder 
approximation. As discussed above, there could also be important corrections to the j3 
function when a is this large. Some evidence that this is the case in the MS scheme 



adopted here can be provided by determining a* using the three-loop MS (3 function [IT 
and setting this value equal to a c . For N = 2, this leads to a new value of Nf G /N w 3.175, 
again shifted less than 20% from the two-loop value 3.929. The shift decreases to about 

15% as N -> oo. 

For Nf < Nf the gap equation gives a dynamical mass E(p), where p is the mag- 
nitude of Euclidean momentum. It has some value E(0) at p = and then falls mono- 
tonically with increasing p. The scale E(0) vanishes continuously in a characteristic 
exponential fashion || as Nf — > Nf G from below (equivalent to a* — > a c from above), 
at fixed A. Since the fermions decouple for p -C S(p), the infrared fixed point is only 
an approximate feature of the theory, useful at momentum scales above the decoupling 
scale. Fortunately, the critical behavior of the theory near the transition is determined 
mainly by the momentum range E(p) <C p ^C A in the gap equation integration, where 
the fixed point is a good approximation. Further discussion of the gluon-induced critical 
behavior may be found in Refs. 0, [13|, [14 . 



Turning now to our study of the role of instantons in the chiral phase transition, 
we will derive a corresponding critical coupling or critical number of flavors Nf , arising 
purely from instanton effects. Comparison with the values from purely gluon exchange 
will then indicate the relative importance of the two effects in the phase transition dy- 
namics. 

A nonzero dynamical mass S(p) in the quark propagator i/[A(p) / ~f fJ- p 11 — E(p)] signals 
chiral symmetry breaking. To determine this two-point function we adopt a formalism 
which self-consistently sums the effects on a fermion propagating through a dilute gas 



of noninteracting instantons, giving a gap equation whose kernel is directly related to the 
single-instanton amplitude. For general N, the gap equation takes the form 



with 

T[a(p)} 



4 e 5/6 e CAf-BAf / - - v 
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(10) 



(N -l)\(N-2)\ 

The numerical factors and a dependence in T[a(p)] arise from the amplitude for an 
instanton of size p, integrated over the other collective coordinates @, |15[ . Gauge field 
and quark fluctuations around the instanton background contribute quantum corrections 
which include logarithms that renormalize the bare coupling to its value at length scale 
p. Depending on the order of the calculation, we end up with the one- or two-loop 
renormalization group solution for a(p), with leftover non- logarithmic terms going into 
the numerical prefactor. 

If the two-loop a(p) is to be used, to make use of the infrared fixed point, then 
the constants B and C should be computed to the same order. They have so far been 
computed only through one loop [0, [H|, where they are B = 0.3595 and C = 2.0706 
in the MS scheme. [] The higher-order fluctuations that generate two-loop running in 
exp (— 2n/a) and one-loop running |16| in (2ir/a) also contribute corrections of 0(a) 
to B and C. We expect that a full two-loop calculation will be scheme independent. In 
the absence of such higher order computations, we will simply take a(p) everywhere in 
the gap equation to be governed by the two-loop (3 function and infrared fixed point. For 
our numerical study, we will use the one-loop, MS values of B and C, assuming that the 
higher order corrections will lead only to 0(1) changes, in particular keeping e B ~ 0(1). 

The function D[pm(p)] contains the mass-dependent factors from the fermionic quan- 
tum fluctuations (left over after the factors containing regulated divergences are absorbed 
into T[a(p)]). Following Ref. || we evaluate the argument of D[pm(p)} using the function 
m(p) derived from E(p) weighted by the fermion zero mode wavefunction, 

m{p) = (Vo|S|^ ) = / dxxf 2 (x/2)E(p = x/p) . (11) 



Here, f(x) is a combination of modified Bessel functions [17] arising from the Fourier 
transform of ipo (x) > 

f{x) = -2I 1 (x)K 1 (x) - 2x(h{x)K (x) - I (x)K 1 (x)) , (12) 
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normalized to /(0) = 1. Its asymptotic behavior is f(x ^> 1) ~ ^x 

For arguments pm(p) ^C 1, the fermion fluctuation factor D[pm(p)] has the expan- 
sion D(x) = x + 0(x 3 \nx). For large arguments, the fermions decouple and D[pm(p)] 

1 In the MS-bar scheme these become B = -0.2917 and C = -1.5114. 



exponentially approaches unity [ 18| , which we must multiply by e B ~ 0(1) to account 
for decoupling the fermion factors in T(a). Lacking a full calculation of the instanton de- 
terminant for massive fermions, we adopt a simple form that interpolates between these 
two limits: 

D{x) = { X n ioTX<eB (13) 

1 1 • e otherwise . 

Spontaneous chiral symmetry breaking corresponds to the existence of nonvanishing, 
energetically preferred solutions to the gap equation (|9|). We will first qualitatively discuss 
the existence of these solutions and then summarize the results of a numerical study. To 
begin, we note that a nonvanishing solution E(p) will have some finite value S(0) at p = 0, 
and then decrease monotonically. This behavior, typical for a dynamical mass in a gauge 
field theory, follows from the structure of the gap equation (||) with the factor f 2 (pp/2) 
decreasing monotonically from unity. The zero mode mass in Eq. (|TTD is approximately 
m(p) ~ E(l/p), since xf 2 (x/2) is peaked around x = 1 and integrates to unity. Thus, 
m(p — > oo) = S(0), while m(p) falls rapidly as p — *■ 0. Although the intrinsic scale A sets 
the solution scale S(0), our interest here is in exploring a possible second order phase 
transition near which the gap equation may dynamically enforce S(0) -C A. 

Assuming this to be the case, the integration over p then breaks naturally into three 
regimes. The ultraviolet regime, < p < 1/A, contributes very little to the integral 
because asymptotic freedom ensures a strong suppression of T[a(p)], while D(pm) Nf also 
remains small. In the intermediate regime, A -1 < p < S(0) _1 , a(p) ranges only from 
roughly 0.78a* to a*, and the infrared fixed point dominates the behavior. The upper 
end of this intermediate regime, p ~ S(0) _1 , should dominate the integral due to the 
polynomial increase of D[pm(p)] while r[a(p)] simply approaches its fixed point value r*. 
In the third regime, as p ^> S(0) _1 the fermions decouple from the fluctuations, ending 
the polynomial increase of D[mp] and the fixed point behavior of a. The third regime 
will affect the critical value of Nf but not the qualitative behavior near criticality. 

We can see this by scaling to dimensionless variables s(y) = £(0) -1 m(?//E(0)) with 
y = E(0)p; note that < s(y) < 1. The gap equation at p = then becomes 

1 - + r, f 1 %[ys{y)] N ^ + rV{y)e B ^4r-, ( 14 ) 

Js(o)/a y 2 Ji y 3 s(y) 

for nonvanishing E(0). Thus it is the intermediate regime that controls the critical 
behavior of S(0)/A. An important feature of Eq. (0) is that since s(y) < 1, it cannot 
be satisfied if T* <C 1, which occurs for Nf close enough to 11N/2 (that is, small a*). 
Only the chirally symmetric solution S(0) = exists for this range of Nf. 

As we decrease Nf to N'j 1 (the critical number of flavors for the instanton kernel), 
T* will eventually reach a critical value large enough to allow nonzero solutions. The 
structure of Eq. ([14]) indicates that this will correspond to maximizing the intermediate 
integral, and therefore to taking its lower limit to zero. For Nf slightly less than Nf , 
E(0)/A must be small, indicating a continuous phase transition. We expect Nf > 3.7N, 



since if a critical value occurs it can only be before T* reaches its maximum as a function 
of Nf, corresponding to a* ~ tt/N. We will not here determine analytically the behavior 
of S(0) as Nf — > Nf . Instead, having qualitatively seen that there exists a continuous 
phase transition at a critical value Nf , we now turn to a numerical study of the transition. 

The numerical results reported in this letter are restricted to the case N = 2. We 
solve Eq. (|)ona one dimensional lattice, iteratively relaxing the discretized S(p) from 
an initial guess to a self-consistent shape, at each stage numerically integrating to get 
m(p). We expect, from our qualitative discussion, that the dominant range of integration 
will be approximately A -1 < p < S(0) _1 . We use the exact solution Eq. (|6|) to the two- 
loop renormalization group equation, for all p < S(0) _1 (more precisely, for p below the 
value at which D(pm) = e B , solved for at each iteration), using the MS value B = 0.3595. 
For p above this fermion-decoupling value, we match a(p) onto the (3 function solution 
for Nf = 0. After a finally grows too large for perturbative running, we simply fix it at 
a constant value a max = 2tt/N; the far end of the infrared range of integration is safely 
subdominant, and we have checked that this approximation is unimportant. 

The result is a shape S(p) for each value of Nf. To study a possible phase transition 
we examine the behavior of E(0), for fixed A, while varying Nf. The numerical results 
confirm the qualitative discussion above: when Nf approaches the critical value Nf « 
4.77N « 9.54 from below, S(0) vanishes continuously in the manner of a second order 
phase transition. Fig. la displays this behavior for A^ = 2; in contrast to the exponential 
behavior arising from gluon exchange, Fig. lb indicates a power law behavior. 

log 10 [S(0)/A] log 10 [E(0)/A] 

2r 




Figure 1 (a): log 10 (E(0)/A) against Nf/N for fixed N = 2. Extrapolation identifies 
the critical value Nf = A.77N. 
(b): Plotting the same data against log 10 (Nf — Nf)/N shows that E(0) 
follows a power law in (Nf — Nf). 

The critical coupling and Nf for an instanton kernel are numerically very similar to 
the values for the gluon-exchange kernel, basically independent of whatever 0(1) values 
are chosen for B and C . In any case, a more complete calculation which combined the two 
kernels would lead to a somewhat smaller combined critical coupling (or larger critical 
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Nf); the qualitative point is that the two effects are rather similar in magnitude. 

Finally, we discuss the validity of the dilute gas approximation, incorporated in the 
gap equation used here to describe the phase transition. It allows the amplitude for a 
fermion propagating in the field of a single instanton to be summed over multi-instanton 
configurations, neglecting instanton interactions. The validity of this approximation de- 
pends on the relative magnitude of the dominant instanton size and the typical separation 
distance between instantons in multi-instanton configurations. As noted above, the dom- 
inant instanton size is of order S(0) _1 , which grows without bound near the critical 
point. But the average instanton separation does just the same, since it is controlled 
by essentially the same instanton amplitude integral, also dominated near p ~ £(0) _1 . 
Crudely estimating the instanton density by 



n 



dp T[a{p)\ 



D[pm{p)} 



Nt 



E(oy 
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we ask that there be fewer than one instanton per instanton four- volume £(0)~ 4 



(15) 
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a * 



2N 



-2-ir/a* 



<1 . 



(16) 



2 this requires a c only slightly smaller than ir/N, that is Nf > 26N/(7—N 



ForN 

3.7 TV. Our numerical result for N = 2 is on the safe side of this limit, giving some 

reassurance that instanton overlap does not violate the dilute gas approximation. For 



larger N, the non-overlap condition (16) becomes 



,5/6 



-N 2 
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2vr 
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,1— n/Na* 



2N 



<1 



(17) 



putting a somewhat stricter bound on a c i and Nf . 

It is clear from this last expression that at large N the validity of the dilute gas 
approximation becomes a delicate matter depending sensitively on the value of a*, and 



also on B and C. If C were neglected in Eq. [17], the remaining expression |Tj| would 
increase exponentially with N, for a* in the range required to trigger chiral symmetry 
breaking. Since the same factor appears in the instanton amplitude entering the gap 
equation, the increase would also affect the dynamics of chiral symmetry breaking. That 
is, the effect of instantons would grow with N, at least up to the point where the dilute gas 
approximation breaks down. Whether this actually happens depends on the prefactor 
term involving B and C . Until these constants are computed through two loops and 
demonstrated to be scheme independent for a(p) = a*, the relevance of instantons to 
chiral symmetry breaking remains uncertain in the large N limit. 

In summary, we have studied the role of instantons in the zero-temperature chiral 
phase transition in SU(N) gauge theories, using the number Nf of fermion flavors as the 
control parameter. The key feature of these theories is that for a range of Nf including 



the critical value for the transition, the two-loop (3 function exhibits an infrared fixed 
point. This allowed us to discuss qualitatively the existence and behavior of solutions 
to a gap equation, whose kernel arose from the propagation of fermions in an instanton 
background in the dilute gas approximation, and to present numerical solutions. We 
found a critical value of Nf, below which chiral symmetry breaking occurs, which is 
comparable to that generated by gluon exchange alone. We conclude that for small N, 
instantons play a role comparable to that of gluon exchange in the chiral phase transition. 
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